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Abstract 


This study investigates the impact of equivalent roughness calculation method on the 
optimal design of composite channel. Seventeen different methods are investigated. Three 
methodologies have been used to solve these models that are division of area into sub areas 
by vertical lies at every places of change in geometry or roughness, division by bisecting 
lines at every places of change in geometry, and incorporating perimeter along the 
imaginary line. Techniques of genetic algorithm and classical method (SQP) have been 
used to solve the optimization problems in all the cases. All these models have been 
solved for three scenarios having no restriction, restriction on top width, and restriction in 
side slopes. In this study, a technique has been proposed to determine the roughness along 
the imaginary vertical lines and this technique has been used in the problem formulation of 
optimal design of trapezoidal composite channel. The results obtained in this study 
indicate that the choice of equivalent roughness formula can impact channel design 
considerably and needs to chose carefully giving due consideration to actual conditions in 
field being closest to the assumption of the method. Also the technique of GA was 
successfully employed for the purpose of optimal design of composite channel. 
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Chapter 1 


Introduction 


1.1 General 


Manmade channels are major conveyance systems for delivering water for various 
purposes such as irrigation, water supply, flood control, inland navigation, and 
rehabilitation, etc. Channel design projects are one of the most expensive civil engineering 
projects in the world. The total cost of these projects includes the cost of construction and 
cost of operation, maintenance and repair (OMR). The lengths of channels carrying water 
for various purposes can be very large (about hundreds of kilometers) so that any savings 
in the channel projects tends to save a large amount of money and manpower. The main 
concern in the design of a channel system is to minimize its construction cost to safely 
carry the design discharge. The total construction cost of channel mainly includes 
excavaition cost and lining cost. Different lining materials for different sides may be 
needed due to one or more of the following reasons: (a) certain site conditions; (b) limited 
right of way; (c) availability of lining materials in local areas. 

Since different materials may be needed in lining of different sides and bed, the resistance 
offered to flowing water along sides and bed in such channels may be distinctly different 
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resulting in different roughness along the wetted perimeter of a charmel cross section. 
These types of channels are called composite channels. In composite channels, the 
roughness and hence shear stress may be distinctly different in different parts of the 
channel cross section, hence an equivalent uniform roughness coefficient is needed. The 
equivalent roughness coefficient equation incorporates flow geometric elements and 
corresponding roughness coefficient value (Chow, 1959). Traditionally, the equivalent 
roughness coefficient of a composite channel cross section is conveniently expressed for 
Manning’s roughness (n) using weighted sum of the entire resistance coefficients along the 
perimeter. In computations, usually finite discretization approach is used by dividing the 
cross-section into number of sub cross-section areas, wetted perimeters, and hydraulic radii 
and the weighing factor is the function of these parameters (Yen, 2002). 

For composite channels, depending upon different assumptions about the relationship of 
discharge, velocities, forces, and shear stresses between component subsections and total 
cross section, numbers of formulae has been proposed. Since each formula has its own 
assumptions and limitations; they are distinctly different from each other. In this study, 
the impact of different assumptions in calculation of equivalent roughness on optimal 
design of composite channel has been studied. In the various formulae, weighing factors 
are calculated by dividing the cross-section into the number of sub-section areas, but the 
location of imaginary boundary line separating different sub areas is not clear. There are 
several assumptions for locations of imaginary boundary line given by Yen (2002). Some 
of these include; (a) vertical lines extending from every break point of the geometry or 
boundary roughness up to the surface water; (b) bisect lines of every angle at the geometry 
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or roughness break point; (c) a horizontal or near horizontal line joining the two breaks of 
the channel at the bank- full stage, separating the compound channel into two parts: the 
lower main channel section and upper flood channel section; (d) a variation of (c) by 
further subdividing the lower main channel and upper main channel by bisect angle lines 
or vertical lines. Since all the previous researchers in the area’ of channel optimization used 
method (a) that is the simplest approach for dividing total area into segments, in this study, 
a detailed examination of the influence of the position of the line of sub area division by 
applying methods (a) and (b) on the optimal design of the composite channel has been 
investigated. 

Some of the formulae are based on the assumption that there are different velocities in 
different subsection. When the velocities in the two adjoining areas are very different, then 
there exists a very high velocity gradient and hence shear stress at the position of 
imaginary vertical line. Therefore, it is needed to modify the segmental wetted perimeter 
to take care of the resistance to flow at such boundaries. This has been investigated in this 
study by incorporating the portion of the imaginary vertical distances corresponding to 
such boundaries in the segmental wetted perimeters. 

The optimization problem for channel design is a c®nstrained optimization problem 
involving minimization of total construction costs. Manning’s equation of uniform flow in 
the form of an equality constraint is normally employed. Any equality constraint that is 
nonlinear in nature would make the optimization problem a nonconvex one (Deb 2001). 
Classical optimization solution procedure may not be suitable to solve such an 
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optimization problem as they have tendency to converge to a local optima. Since all the 
classical and gradient based search techniques start search from a single initial solution, 
there is higher probability to getting stuck into a local optima. Hence, this type of problem 
must be solved with a global optimization strategy that is based upon: (a) combination of 
deterministic and probabilistic approaches; and (b) systematic evolution of the solution in 
the direction of the global improvement, and concept of competitive evolution (Duan et al. 
1993). Therefore, evolutionary method of optimization (Genetic Algorithm) is used to 
solve the problem as it uses probabilistic and deterministic approach starting its search 
from a set of possible solutions and has a higher probability of obtaining global optimal 
solution. 

1.2 Objectives of the thesis 

The objectives of this study are as follows: 

1. Analyzing the impact of equivalent roughness coefficient method on the optimal 
design of composite channel 

2. Influence of the position of sub area division on the optimal design of the channel 

3. Developing a new optimization program based on the incorporation of the partial 
wetted perimeter along the imaginary vertical boundaries as area subdivision lines 

4. Evaluating' various scenarios of optimal design based on the practical situation. 
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1.3 Organization of the thesis 


This chapter presents the brief overview of the problem of the optimal design of composite 
channels. Need of using genetic algorithm to solve the optimization problem of composite 
channel cross-section design has been also illustrated here. This chapter also presents the 
need of different scenarios in optimal design of composite channel, influence of line of 
sub-area division and also accounting the effect of wetted perimeter along the imaginary 
subdivision on optimal design cost. Chapter 2 consists of literature review of the optimal 
design of open channel cross section and use of genetic algorithm in past in the closely 
related areas. Chapter 3 gives the brief overview of binary coded genetic algorithm that 
has been used in this work. Chapter 4 describes different models, scenarios, approaches, 
and the solution methodology to solve the optimization problem of composite channel 
cross section design. Chapter 5 presents the results and their discussion and finally 
summary, conclusion, and scope for future work have been discussed in chapter 6, 
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Chapter 2 


Literature review 

2.1 General 


Channel design projects have been in existence since earlier civilization. The time, when 
man started agriculture, which is rain dependent and rainfall is highly uncertain, one 
needed artificial methods to irrigate crops. That time artificial channel designing projects 
came into existence. There are a lot of evidences that in ancient river valley civilizations 
(Sindhu River Valley Civilization, 3500-2500 BC) man had designed artificial channels. 
In ancient literatures like Arthashastra by Kautlilya (about 200 BC) one can find methods 
to design channels. But in modem era, the scenario has been completely changed. 
Manmade channels are used not only serving for the irrigation purposes; but also for fresh 
water supply, inland navigation, flood control rehabilitation etc. With the development of 
sciences and optimization techniques researchers and hydraulic engineers have shown 
interest in optimizing the channel design problems. Initially, methods of calculus were 
used to obtain most efficient hydraulic section with least resistance. 
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2.2 Classical methods 


In last 25 years, researchers have done a lot of wok in the area of optimal channel design. 
Mironenko et al. (1984) proposed that parabolic shaped unlined and lined canals are 
hydraulically more stable. The extra stability of the parabolic channels is due to the fact 
that the only one point of the cross section is at the maximum side slope. A lined parabolic 
channel has no sharp edges of stress concentration where cracks may occur. Guo and 
Hughes (1984) designed optimal channel cross-section from the first principle of calculus. 
They presented an analytical procedure to determine the dimensions of a trapezoidal open 
channel. They proposed two strategies; the first strategy yields optimal chaimel shape that 
provides the least hydraulic resistance and the second one yield least cost channel 
geometry. Differing from tradition, for the first time they incorporated freeboard as input 
parameter in optimization problem, as in traditional methods, freeboard was provided after 
getting optimal dimensions. They presented the results in a graphical form for least 
hydraulic resistance using a series of figures that can be used in design by practicing 
engineers. The results revert back to traditional solutions if freeboard parameter is set 
equal to zero. They did not present any graphical solution for minimum construction cost 
but proposed a cost function, related to excavation depth and channel width that yields a 
general expression for a least cost trapezoidal cross-section. Solutions for least cost 
rectangular channel cross-section with or without freeboard were also demonstrated in 
their study. 
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Flynn and Marino (1987) presented the concept of economically optimal section. They 
constructed canal cross-section designs for uniform flow and compared by non- 
dimensional shape parameter. They exposed basic relationship among the cross-section 
shapes and design variables that are wetted perimeter, water depth, water surface width, 
cross section area, lining volume, and the excavation volume. The relations are used to 
uncover robust rules that can determine optimal cross-section of canal design for 
elementary problems directly from flow information such as capacity, velocity, slope and 
roughness. For problems involving complex limits and economics, the relations are 
combined with optimization methods to solve for the economically optimal cross section. 
Loganathan (1991) presented optimality conditions for a parabolic channel cross-section. 
He accounted for freeboard and put limits on velocity and channel dimensions in his 
optimization formulation. He stated analytically the necessary conditions for minima. He 
found that formulation with an objective function consisting of minimization of excavation 
and lining costs might favor slightly wider channels when the ratio of cost of lining to cost 
of excavation is high. He has given design aids in the form of tables along with formal 
procedures to create additional tables, which are necessary for dealing with depth 
dependent freeboard. 

Froehlich (1994) implemented, for the first time, limited top width and depth as additional 
constraints in cross-section of nonerodible trapezoidal open charmel optimization problem 
that are likely to occur in the field. He solved the constrained minimization problem of 
finding the best hydraulic channel cross-section using Lagrange multiplier method with the 
dimensionless forms of the constraints and objective fimction. Monadjemi (1994) used 
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Lagrange’s method of undetermined multipliers to determine the best hydraulic channel 
section for a given flow roughness and bed slope. He solved the optimization problem of 
the channel section by minimizing the wetted perimeter or cross-section subjected to a 
constraint, which may be any flow equation (i.e. Manning’s equation of uniform flow). He 
applied this method on the standard sections as well as the roimd bottom triangular 
sections and found that, the minimizations of cross-section area as well as minimization of 
wetted perimeter are mathematically equivalent. Swamee (1995) investigated explicit 
equations for the design variables of various canal sections. He found that triangular and 
rectangular sections are equally optimal in the sense that they yield the same minimum 
area and semicircular section has least flow area and flow perimeter. Thus, constrained by 
straight line geometry, the optimal canal section becomes closest to the semi-circular 
section, Swamee et al. (2000) investigated the optimization problem of open channel 
cross-section design considering minimum seepage loss. They applied nonlinear 
optimization technique to solve the explicit equations for design of minimum seepage loss 
channel section that are obtained by using seepage loss equation and uniform flow 
equation. 

As lining materials is one of the basic constituents in the channel design cost, a lot of 
attention is needed to select the lining material for the whole wetted perimeter. There may 
be several site conditions, for that it becomes necessary to use different lining materials on 
different sides and bed. Cost of lining materials mainly depends upon their availability in 
local area otherwise transportation cost will increase the whole construction cost. 
Attention on this phase of the problem was given by Trout (1982) who developed direct 
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algebraic techniques to determine open channel cross-section design which minimizes 
lining material cost when different materials are used for lining bed and sides. He 
presented solutions in the form of graphs, which indicates optimal parameter combination, 
and cost of deviation from the optimal design. He considered only minimization of the 
lining cost of the trapezoidal channel but did not consider minimization of the cross- 
sectional area. Das (2000) formulated optimal design problem in nonlinear optimization 
framework with cost per unit length of composite channel as objective function. 
Manning’s equation, positive values of design variables and specified values of side slopes 
or top width were used as constraint. He considered the excavation cost and lining cost as 
objective function. He used Lagrange Multiplier method to convert the constraint 
optimization problem to unconstraint optimization problem and developed a computational 
methodology using first order necessary condition for optima to solve this problem. 
Recently, Jain et al (2004) solved the problem of composite channel design using genetic 
algorithm for first time for different scenarios including top width constraint, side slope 
constraint and velocity constraint. 

There are many equations that can be used to calculate the equivalent roughness of 
composite channels. Yen (2002) presented a list of many of these equations that can be 
used to formulate different models for composite channel optimization problem. There are 
several inherent assumptions and limitations to each equation of equivalent roughness that 
are also applied to optimization models formulated using these equations. The final results 
may be different for different models because of these different equations and their 
assumptions. In the previous approaches. Das (2000) used Horton’s (1933) equation for 
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calculation of equivalent roughness which is a lumped approach and considers equal 
velocity in all segmental areas of the composite channel that is equal to mean velocity of 
the channel. Some other methods of calculation of equivalent roughness also assume same 
approach in the sense that there is a constant velocity across the cross-section but few of 
them allowed the spatial variation across the cross section. It should be realized that the 
different roughness coefficient across the channel cross section should modify the velocity 
distribution. This drawback in the Horton’s method can be overcome by using Letter’s 
(1933) approach of calculating equivalent roughness in composite channels. This method 
does not assume segmental velocities as equal to mean velocity and allows spatial 
variations in segmental velocities across the cross section. Jain et al (2004) proposed a 
new model that was used in optimization problem considering the spatial variation across 
the cross section. 

2.3 Genetic Algorithm 


In last decade, the researchers working in the area of hydrology and water resources have 
shown interest in using GA to solve optimization problem. Wang (1991) successfully 
applied GA in calibration of conceptual rainfall runoff model for data from a particular 
catchment. There were seven parameters in his model and all seven parameters have been 
optimized by using GA. Cieniaveski et al. (1995) used GA to solve multi-objective ground 
water monitoring problem and compared the results with the results from simulated 
annealing. They found that GA is able to find a larger number of non dominated points on 
the trade off curve while simulated annealing found fewer points. Liang et al. (1995) used 
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GA in search for optimal values of catchment’s calibration parameters. They linked GA to 
a widely used catchment model, the storm water management model (SWMM) and applied 
to a catchment in Singapore. They considered six storms and out of them three were used 
for calibration and three were used for verification of the results. Their study showed that 
GA required only a small number of catchment model simulations and yet yielded 
relatively high peak flow prediction accuracy. Gentry et al. (2001) employed GA to 
determine areas contaminating a semi confined aquifer. There are many other researchers 
who used GA successfully in different areas of water resources such as urban drainage 
modeling (Rouch et al., 1999); optimal design and rehabilitation of water distribution 
system (Wu, 2001); estimation of aquifer parameters (Samuel, 2003); location of aquifer 
leakage areas (Gantry, 2003). Jain et al. (2004b) used real coded GA to train ANN rainfall 
runoff model and Jain et al. (2005) used real coded GA to determine unit pulse response 
fiinction using historical data from watershed. There are many other works that have been 
carried out by different researchers showing that GA can be efficiently used in the complex 
problems of water resources. Optimal channel designing problems are not far away from 
the reach of GAs. Recently, Jain et al (2004) solved the problem composite open channel 
design using genetic algorithm for first time for different scenarios including top width 
constraint, side slope constraint and velocity constraint. They also proposed a model based 
on the distributed approach, which allows spatial variation in velocity across a composite 
channel. They presented results in form of optimal design graphs to determine dimensions 
of optimal trapezoidal composite channel section. 
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2.4 Summary 


Channel designing projects started in very old ages but the main concern in this area was 
given in last two decades. Initially, optimal design was carried out with the methods of 
calculus to obtain most efficient section with least resistance. The techniques employed to 
solve the problem of optimal channel design mainly has been classical method; however, a 
few researches have used GAs for this purpose. The channel optimization problem 
involves minimization of total construction cost with uniform flow equations and non 
negativity constraints. Various scenarios involving additional constraints have been 
attempted. Recently, some researchers have proposed solutions to the optimal design of 
composite channels. This study extends the existing work in the area of optimal design of 
composite channels. 
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Chapter 3 

Genetic Algorithm 


3.1 Need of Genetic Algorithm 


The need for solving optimization problems arises in almost every field and, in particular, 
is a dominant problem in water resources systems. Consequently, many analytical and 
numerical optimization techniques have been developed. There exist a great number of 
functions, that are discontinuous, nondifferentiable, nonconvex, and multimodal, which are 
beyond analytical methods and pose profound difficulties for traditional numerical 
techniques. Generally, traditional optimization problem solution techniques are based on a 
deterministic relationship between the model parameters and its performance but now a 
days these techniques have been unable to optimize such complex problems. (McKinney 
and Lin 1994, Reed et al. 2000). In addition to these, there are some common difficulties 
with most of the classical direct and gradient based techniques; (a) convergence to an 
optimal solution depends on the chosen initial solution; (b) most algorithms tend to get 
stuck in a sub optimal solution; (c) an algorithm efficient in solving one optimization 
problem may not be efficient in solving a different optimization problem; (d) algorithms 
are not efficient in handling problems having a discrete search space; and (e) if search 
space contains some local optima, classical methods get attracted to any of locally optimal 
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solution, and there is no escape. Therefore, new and more robust optimization techniques, 
which are capable of handling such problems, are needed. The search for more efficient 
numerical optimization methods has led researchers to reproduce some mechanisms that 
are naturally robust such as genetic algorithms (GAs). 

GA is a search technique based on the concept of natural selection inherent in the natural 
genetics, which combines an artificial survival of the fittest with genetic operators 
abstracted from nature (Holland 1975). The major difference between GA and other 
gradient based or traditional optimization solution techniques is that GA start search with a 
population of many possible solutions; whereas other classical optimization solution 
techniques start their search with a single point. Hence, GAs are less susceptible to getting 
“stuck” at local optima than gradient search or other traditional classical optimization 
methods, rather they search for global optima. In the past decade, GAs have been applied 
successfully to a number of problems such as optimizing simulation models, fitting 
nonlinear curves to data, solving systems of nonlinear equations, and machine learning 
(Deb, 2001). 

3.2 Basics of Genetic Algorithm 


First of all, an initial population of solutions is generated using a random generator. Every 
population is presented by a set of parameter values that describe the problem. 
Traditionally, GAs have been developed by using binary coding, in which a chromosome 
or string, is represented by a string of binary bits i.e., O’s and I’s which are called genes 
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that can encode integers, real numbers, or anything else appropriate to the problem. Binary 
strings are easy to operate on, and within any gene, binary representations can be mapped 
to values in the range feasible for the variable represented (Goldberg, 1989). In this study, 
the parameters or decision variables are encoded as substrings of binary digits having a 
specific length. The length of the substring is determined according to the desired solution 
accuracy, and is dependent on the specified range of a parameter i.e., upper and lower 
bounds and the precision requirement of the parameter. By increasing the substring 
length, more accuracy in the solution can be achieved. Any precision for variables can be 
achieved by changing the string length and lower and upper bounds. 


z, = x,™'" + 


max __ min 


2 ‘‘ -1 


DV(s,) 


(3.1) 


Where It is length used to code the z-th variable and DV(Si) is the decoded value of the 
string Si, xi is actual value of the variable, Xi """ and Xi are lower and upper bounds for 

the variable and the complete string will be 5 = Y." ^ Sj . The above mapping function 


represented by equation (3.1) allows any arbitrary precision for a variable and different 
precision for different variables, by changing the length of substring. 


3.3 Fundamental Operations of GA 


GA starts with randomly generating an initial population of possible solutions. For most 
GAs these candidate solutions are represented by chromosomes coded using a binary 
number system. The GA that employs binary strings as its chromosomes is called binary 
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coded GA. These chromosomes are evaluated based upon their fitness function, in terms 
of certain objective function. The binary coded GA consist of three basic operators, 
selection operator, crossover operator and mutation operator. These are discussed in the 
following sections. 

3.3.1 Fitness Function 

In GA, solution to the problem is evaluated using fitness functions, also known as 
objective function. The evaluation of function determines the cost of the solution, which 
includes different costs involved in the objective function. Based upon the results of such 
functions, evolutionary procedure may be applied to a set of solutions. In case objective of 
optimization is maximization of the objective function, a solution with larger fitness value 
compared to other solutions is better but if the objective of optimization is minimization of 
the objective function, a solution with smaller fitness value is better as compared to other 
solutions. 

3.3.2 Selection Operator 

The objective of selection operator is to make duplicates of good solutions and eliminate 
bad solutions in the population while keeping the population size constant. There are a 
number of methods Uke tournament selection, proportionate selection, and ranking 
selection by which this task can be achieved, which include identification of good 
solutions in population, production of multiple copies of good solutions, and elimination of 
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bad solutions from the population so that bad solutions can be replaced by multiple copies 
of good solutions. These can be found in Deb (2001). 

3.3.3 Crossover Operator 

Crossover and mutation are driving force in any evolutionary algorithm. The theory 
behind the crossover operation is that by exchanging important building blocks between 
two strings that perform well, the GA attempts to create new strings that preserve the best 
material from two parent strings. Building blocks are the chromosomes that contain 
information in the form of genes. The number of strings in which material is exchanged is 
controlled by crossover probabiUty. Crossover operator is a recombination operator that 
picks up two chromosomes, called parent strings, in mating pool. Crossover location is 
selected and the genetic information is exchanged between the parent chromosomes with a 
certain probability (pc) to create two new chromosomes, called child strings. Crossover 
probability can lie between (0-1.0) and is generally taken from range (0.6 -1.0). This 
procedure is repeated to other pairs of chromosomes to give a new generation of child 
population having better fitness value than that of the parents. Crossover has a wide range 
of possible types, viz., one-point, multiple-point, uniform, intermediate arithmetical, and 
extended arithmetical which can be found in Deb (2001). The functioning of a single 
point crossover is explained in figure (3.1) in which Pi and P 2 represent the parent 
chromosome and Ci and C 2 represent the two children. The Character represents the 
location of single point crossover in both parents. 
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Figure 3.1: Crossover operator 


3.3.4 Mutation Operator 

GA is a very aggressive search technique and discards the chromosomes having poor 
fitness values which may cause our search to converge to local optima. To eliminate this 
possibility and get global optima, mutation operator is used. Mutation provides a small 
amount of random search, and helps to ensure that no point in the search space has a zero 
probability of being examined. The mutation operator generally facilitates the child 
population to converge towards optimal solution even if the initial population is far from 
optima. The bit wise mutation operator in binary coded string changes a 1 to a 0 and vice 
versa with a small mutation probability of pm. It is used to maintain the diversity among 
the child population from one generation to another generation. Proper selection of 
mutation probability is important for proper working of GA. In previous studies, it has 
been suggested that a mutation probability inversely proportional to the population size is 
enough to prevent the search to converge into a local optima. 







20 


Figure (3.2) describes the process involved in mutation. “P” represent the parent 
chromosome and “C” represent the child chromosome. Character represent the 
location of mutation where bit 1 changes to 0. Mutation probability can lie between (0- 
1 .0) but generally taken very small. 



11 0{100101 


O 


1 1 0{000 1 0 1 



Figure 3.2: Mutation operator 


The process of selection, crossover and mutation is repeated from one generation to the 
next, achieving better fitness in each generation. The search is terminated using certain 
stopping criteria in terms of the convergences of the solution. Figure (3.3) shows the 
flowchart for a simple GA. 
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Chapter 4 

Model Development 


4.1 Introduction 


In this study seventeen nonlinear optimization problems for optimal design of trapezoidal 
composite channel cross section have been investigated. All these problems are 
constrained minimization problem and consist of same objective function which is cost of 
construction of the trapezoidal composite chaimel. Uniform flow conditions are ensured 
using Manning’s equation as a constraint. The total cost of construction includes cost of 
excavation that depends upon cross sectional area, and cost of lining that depends upon the 
whole wetted perimeter. In this study, attempt has been made to design composite channel 
in which different lining materials have been used on different sides and bed and having 
different costs per unit length. Basic differences in all these optimization problems are that 
they are subjected to different approaches of calculation of equivalent roughness. Since 
roughness is changing along the wetted perimeter and is distinctly different in different 
parts; an equivalent roughness («c) is required to be used in Manning’s uniform flow 
formula (Chow 1959). There are various methods of calculation of equivalent roughness 
that have been developed by various researchers in the past. The basic differences in these 
methods are that they have been developed based upon certain assumptions and all 
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methods have their own limitations. When these equations are used in optimization 
problem, their assumptions and limitations are liable to modify the entire search space of 
the problem. This study investigates the impact of roughness on the optimal design of 
composite channel. The basic formulations of all these optimization problems are as 
follows: 


4.2 Problem Formulation 


The objective function is to minimize cost function, which is cost of construction per unit 
length of channel, consist of excavation cost and lining cost and subjected to uniform flow 
condition and non-negative parameters of design. Before discussing the problem 
formulation in detail let us define the properties of trapezoidal chan n el cross section that 
will be used in problem formulation. 

4.2.1 Properties of Composite Trapezoidal Channel 

Figure 4.1 shows a composite trapezoidal channel cross section. It has base width b, depth 
of water y, freeboard / and side slopes zy.l (Horizontal: Vertical) and zy.l (Horizontal: 
Vertical). Let Ty, be the flow top width and 2/be the total top width of the cross section. 
The case when vertical lines are imaginary lines that divide the total area into sub areas as 
in figure 4.1 let Awu Av ,2 and Ay,} be the wetted flow subsection areas and Ay, be the total 
wetted flow area. Similarly, let P^be the flow wetted perimeter of the cross section and 
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P-wh Pw2 snd Pw3 be the wetted lengths on each side. The hydraulic parameters can be 
calculated using the following equations: 



Figure 4.1: Trapezoidal composite channel section 




(4.1) 

(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 

(4.7) 

(4.8) 
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Similarly, subsection areas Ai, A2, and A3, the total area At, total wetted perimeter Pt and 
the total channel cross section top width Tt of the trapezoidal channel cross section 
accounting for freeboard can be written as follows: 


A =:EAy±iy. 

‘ 2 

A =hizil)L 
' 2 

^3=b(y+f) 

A,=b{y+f)+{z^ +Z2) 

P,=y^l 

Tt=b + {z^+z^){y+f) 


(4.9) 

(4.10) 

(4.11) 

(4.12) 

(4.13) 

(4.14) 


Let, Pi and P2 are the length of sides including the freeboard and P3 be the wetted bed 
width of the channel. These can be written as: 

p,=[V?^)(y+/) (4.15) 

P,=[Vz,^+lj(y+/) (4.16) 

P,=b 


(4.17) 
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But in case when the bisecting lines are dividing the total area into sub areas as in figure 
4.2, subsection areas A^,j, and^f^j becomes modified and their equations will be 




(4.18) 


-\ 4^1 


A.3 =by + {z, + 22 )^ — 


(4.19) , 

(4.20) 


Method of calculation of these areas is given in Appendix A. 


4.2.2 Uniform Flow Equation 


The Manning’s equation of uniform flow in a channel can be expressed as: 


4^0 P"" 


(4.21) 


Where Q is the total discharge in the channel (mVs), So is the bed slope of the channel, n is 
the roughness coefficient of the channel cross section, A is the channel cross sectional area 
(m^) and P is the total wetted perimeter (m). 
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4.2.3 Basic Problem Formulation 


The optimization problem can be represented by the following equations: 


Minimize: 


C(b, z,,Z2,y) = c^A,+c^P^+c^P^+c,P^ 


(4.22) 


Subjected to: 


0(Z>, y, z,, Z 2 ) = 


4 ^ 


(4.23) 


b> 0 ; z/ > 0; > 0; y > 0 


Where At is the total cross sectional area that is area of excavation as in equation (4.12), 
Pi, P 2 and P 3 are length of sides including freeboard and defined in equations (4.15), (4.16) 
and (4.17), respectively, c; is cost of excavation per sq. meter, C 2 is cost of lining per meter 
length of perimeter Pi , cs is cost of lining per meter length of perimeter P 2 , C 4 is cost of 
lining per meter length of perimeter P 3 ,ni, n 2 and «5 are Manning’s coefficients for 
different lining materials used on different sides of perimeter and «c is equivalent 
roughness. 

Since the flow equation makes an equality constraint and while solving with GA it 
becomes a problem because GA has its inherent limitation of not being able to handle 
equality constraints; therefore it becomes necessary to convert the equality constraint into 
inequality constraint. It was achieved by introducing a dummy variable e in the 
optimization problem formulation. Hence, the constraint function will be 
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0(6, y, Zp Zj) = 




>0 


(4.24) 


4.3 Equations and Assumptions of Different Models 


Seventeen models based upon methods of calculation of equivalent roughness have been 
formulated. Basic formulations of all these optimization models are as given above. 
Objective functions of all these models are exactly same (equation 4.22) but the basic 
difference is methods of calculation of equivalent roughness that makes modifications in 
search space. Methods of calculation of equivalent roughness and their inbuilt 
assumptions and limitations are as follows: 


Model A 

I 


This method of calculation of equivalent roughness was given by US army corps of 
engineers (USCOE), Los Angles (1973). In this method equivalent roughness is calculated 
by summing component n weighted by area ratio. The shear assumption involved in this 
method is that the shear velocity is weighted sum of sub area velocity i.e 




fv.] 



[vj 


UJ 


(4.25) 


(4.26) 


Equivalent roughness in this method will be as follows: 
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(4.27) 


Model B 


In this method, equivalent roughness is calculated by taking square root of sum of the 
square of the component n weighted by sub area ratio. Assumption involve in this model 
is total resistance force is equal to sum of sub area resistance forces i.e. 


PrRS = YuP<yR<Si 


and 


fK] 



2/3 


(4.28) 

(4.29) 


Equivalent roughness in this method will be as follows: 



Model C 


(4.30) 


In this method equivalent roughness is calculated by dividing total area by sum of sub 
areas normalized by component equivalent roughness. Assumptions involved in this 
method are that total discharge is sum of sub area discharges i.e. 

Q = = = (4.31) 
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and 




(r] 





(4.32) 


Method of calculation of equivalent roughness in this method will be as follows: 


A 


(4.33) 


Model D 


This method of calculation of equivalent roughness was given by Colebatch (1941). In this 
method also sub area ratio is used as weighing factor for component n to calculate the 
equivalent roughness. Assmnptions in this method are that the total cross sectional mean 
velocity is equal to sub area velocity i.e. 

F = F, (4.34) 

and energy grade line for total cross sectional area is parallel to that of sub section areas 
i.e. 


S = S, 


(4.35) 


Method of calculation of equivalent roughness in this method will be 


”, = 


. 2/3 


(4.36) 
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Model E 


This method of calculation of equivalent roughness was given by Horton (1933). In the 
above four methods (A to D) component n were weighted by sub area ratios to calculate 
the equivalent roughness but in this method sub area wetted perimeter ratio is used as 
weighing factor to calculate the equivalent roughness. Assumptions in this method are 
same as above in model D as in equations (4.34) and (4.35). Method of calculation of 
equivalent roughness in this method will be as: 


n 


c 





n2/3 

y 


(4.37) 


Model F 


This method of calculation of equivalent roughness was given by FeUcel (1960). In this 
method of calculation of equivalent roughness the total flow area wetted perimeter is 
divided by sub area wetted perimeters normalized by component n. Assumptions of this 
method are that the total discharge is some of sub area discharges as in equation (4.31) and 

10/3 

(4.38) 





UJ 




Method of calculation of equivalent roughness in this method will be as follows: 
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(4.39) 


Model G 


This method was given by Pavlovaskii (1931). In this method equivalent roughness is 
calculated by taking square root of siim of the square of the component n weighted by sub 
area wetted perimeter ratios. Assumption involve in this model are as follows: 

Total resistance force is equal to sum of sub area resistance forces as in equation (4.38) 

and = (4.40) 


Method of calculation of equivalent roughness in this method is 



(4.41) 


Model H 


This method of calculation was given by Yen (1991). In this method equivalent roughness 
the calculated by summing the component n weighted by sub area wetted perimeter ratios. 
Assumptions involved in this method are that the total shear velocity is wetted sum of sub 
area shear velocities or contributing component roughness is linearly proportional to 
wetted perimeter as in equation (4.25) and the velocity ratio is related to hydraulic radii 
with the equation (4.40). In this method equivalent roughness will be 
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= 




(4.42) 


In the models A to model D sub area ratios are used as weighing factor and from models E 
to H sub area wetted perimeter ratio is used as weighing factor for component n to 
calculate the equivalent roughness. From model I onwards sub area ratio and wetted 
perimeter ratio both are used as weighing factor. Here, another parameter called hydraulic 
radius that is area divided by wetted perimeter is also used. 


Model I 


In this method of calculation of equivalent roughness square root of suimnation of square 
of component n weighted by sub area wetted perimeter ratio and also hydraulic radii ratios 
is taken. Assumptions involved in this method are that the Total resistance force, F, is 
some of sub area resistance forces, ^ F, as in equation (4.28) and total cross sectional 

velocity is equal to sub area velocities as in equation (4.34). In this method equivalent 
roughness will be 


«c = 




^ nfP, 


. 1/2 


(4.43) 
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Model J 


In this method also equivalent roughness is calculated by taking square root of summation 
of squares of component n weighted by sub area wetted perimeter ratios and hydraulic 
radius ratios in a different way. Assumptions of this method are that the total resistance 
force is some of sub area resistance forces as in equation (4.28) 

and {V,IV) = {RJR)"^ (4.44) 


Method of calculation of equivalent roughness is as follows: 




S ./2 


PR 


2/3 


(4.45) 


Model K 

This is another method of calculation of equivalent roughness. This method allows spatial 
variation in velocity across the cross section. It has assumptions that the total discharge is 
equal to sum of sub area discharges as in equation (4.31) and also sub area grade lines are 
inversely proportional hydraulic radius 

= (4.46) 
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This method was given by Letter (1933). This method also allows the spatial variation of 
velocity across the total cross section. Assumptions of this method are that the total 
discharge is equal to sum of sub area discharges as in equation (4.31) and energy grade 
lines for sub areas and total cross sectional areas are parallel as in equation (4.35). 
Equivalent roughness in this method will be 



(4.48) 


Model M 


This method was given by Ida (1960). In this method of calculation of equivalent 
roughness the assumptions are, the total discharge is equal to sum of sub area discharges as 
in equation (4.3 1), and energy grade lines for sub areas and total cross sectional areas are 
parallel as in equation (4.35) and also 
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(4.49) 


Equivalent roughness in this method will be 




Z 


p>R; 


5/3 


(4.50) 


Model N 

This method was given by Yen (1991). This method does not allow the spatial variation in 
the velocity across the cross section. The assumptions involved in this method of 
calculation of equivalent roughness are that the total shear velocity is weighted sum of sub 
area shear velocities as in equation (4.25) and velocity in total cross sectional area is equal 
to velocity in sub areas as in equation (4.34). Equivalent roughness in this method will be 
as: 


= 




P/R 


1/6 


(4.51) 


Model O 


This is another method of calculation of equivalent roughness. Assumptions involved in 
this method are, total shear velocity is equal to sum of sub area shear velocity as in 
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equation (4.25) and velocity is related to hydraulic radius as in equation (4.29). Equivalent 
roughness in this method will be 


= 




PR 


1/2 


(4.52) 


Model P 


This method of calculation of equivalent roughness was given by Yen (1991). This 
method of calculation of equivalent roughness is based upon some assumptions that are 
total shear velocity is equal to some of sub area shear velocities as in equation (4.25) and 
velocity is related to hydraulic radius as in equation (4.44). Equivalent roughness in this 
method will be: 


p^l/3 


1/3 


(4.53) 


4.3.17 Model Q 


Model Q in this study is the one developed by Jain et al. (2004). They proposed a new 
model formulation for optimal design of composite chaimel. They used the same objective 
function that is used in this study but proposed a new constraint which is based upon the 
spatial variation of velocity across the cross section and used the bisection method for area 
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sub division. The constraint function for the optimal design of composite channel 
employed was; 



(4.54) 


All of the above problem formulations are called scenario-I in this study in which there are 
no additional constraints. Some times, it becomes necessary to impose additional 
constraints in the optimization problem to simulate certain site conditions. Two additional 
scenarios were investigated for all of the above models. These scenarios as follows: 


4.4 Scenario 11 

Sometimes for practicing engineers, it becomes a problem that there, is limited width 
available to construct the channel. In order to simulate such type of conditions when the 
top width is limited and restricted up to a maximum value Tmax, as used by Froehlich 
(1994), Das (2000) and Jain et al. (2004), an additional constrain can be imposed on the 
optimization problem. The additionaT constraint to the optimization problem will be as 
follows: 


*^2(^»T>^1>'^2) ~ ^max 7)^0 


(4.55) 
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Where <I>2 is an additional inequality constraint function limiting the total top width of 
composite channel to Tmax. 


4.5 Scenario III 


Sometimes side slopes becomes a problem for practicing engineers because all soils have 
some stability criteria that should be satisfied. For the case when it becomes necessary to 
restrict side slope up to a minimum value and a maximum value due to certain reasons as 
slope stability criteria, or limited right of way, etc, additional constraints can be imposed 
on the optimization problem. The additional constraints to the optimization problem will 
be as follows: 


^3(^1)“^! ^min ^ ® 

(4.56) 

<1)4(22) = Z 2 -Z„i„ ^0 

(4.57) 

*^5(^1) “ ^max ~^i ^ 0 

(4.58) 

‘I>6(22)=-Zmax-^2 ^0 

(4.59) 


4.6 Solution Methodology 

In this study seventeen different models (model A to model Q) have been formulated for 
optimal design of trapezoidal composite channel. Here, all these models differ from each 
other depending upon the method of calculation of equivalent roughness. All the problems 
and scenarios have been solved using two different methods: (a) binary coded genetic 
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algorithm (GA); aird (b) sequential quadratic programming (SQP) from MATLAB-7.01 
optimization toolbox. Results are presented in form of tables in the next chapter. All the 
models have three sceiiarios, namely no restriction (scenario-I), restriction on top width 
(scenario-II), and restriction on side slopes (Scenario-Ill). The values of various defined 
parameters in this study are presented in Table (4.1). These values of various parameters 
were chosen for comparison purposes because these have been used by the other 
researchers in the past. 


Table 4.1 Different parameters used to formulate optimization problem 


Discharge 

(0m"/s 

Freeboard 

(meter) 

Slope 

Cl 

C2 

C3 

C4 

ni 

n2 


(m) 

100 

0.5 

0.0016 

0.6 

0.2 

0.25 

0.3 




H 


The flowchart of working of GA has been presented in Figure 3.3. While using GA to 
solve these problems, we require some user defined parameters. The values of different 
parameters used to solve these problems are listed in Table 4.2. 


Table 4.2 Different parameters used to solve optimization problem with GA. 


Population 

size 

Crossover 

probability 

Mutation 

probability 

Bit length 

used for 

different 

variables 

Total bit 

length 

Number 

of runs 

200 

0.876 

0.0125 

20 

80 

20 
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4.61 First Methodology: Sub area division by vertical lines 


Almost all the researches of the past relating to optimization problem of trapezoidal open 
channel cross section are based upon formulation of problem by the simplest approach of 
division of total cross sectional area into sub areas, that is vertical lines extending from 
every break point of the geometry or boundary roughness up to the water surface as shown 
in Figure 4. 1 . This methodology has been applied to all the models (model A to model Q) 
for ail three scenarios and results by genetic algorithm are presented in the Tables 5.1 to 
5.3 for Scenario-I, Scenario-II Scenario-Ill respectively. Results by SQP (MATLAB 7.01) 
are presented in Tables 5.4 to 5.6 for Scenario-I, Scenario-II, and scenario-III respectively. 

4.6.2 Second Methodology: sub area division by bisecting lines 

As almost in all the models (excluding models E to model H) equivalent roughness 
depends upon the sub areas to total area ratios, method of division of total area into sub 
areas modifies this ratio and hence the equivalent roughness coefficient. Since the method 
of division of total area into sub areas modifies the equivalent roughness coefficient «c, it 
causes modification in the search space of the optimization problem and hence the 
optimization problem will be different. In this study, another approach that is dividing 
total area into sub areas by lines bisecting every angle at the places of changes in geometry 
has been applied to all the models for optimal design of composite channel cross section. 
This method is shown in Figure 4.2. 
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Figure 4.2: Trapezoidal channel section with bisecting lines of angle 

The results of this methodology by Genetic Algorithm are presented in Tables 5.7 to 5.9 
for Scenario-I, Scenario-II and Scenario-Ill respectively and the results by SQP (MATLAB 
7.01) are given in Tables 5.10 to Tables 5.11 for Scenario-I, Scenario-II and Scenario-Ill 
respectively. 

4.7 Proposed Optimization Model 

The division of trapezoidal cross-section into sub-areas can be done in many ways as 
explained earlier. The simplest approach is to use vertical lines at the intersection of side 
slopes and/or change in roughness in the bed. However, when the velocities in the two 
adjoining areas are very different (as is the case looking at the results), then there exists a 
high velocity gradient at the assumed vertical imaginary boundary. The high velocity 
gradients at the imaginar y vertical boundaries call for their inclusion into the total wetted 
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peiimeter. This is because the high velocity gradients at the imaginary vertical boundaries 
essentially lead to an existence of shear stresses at such boundaries. Therefore, we need to 
modify the segmental wetted perimeters in different segments to take care of the resistance 
to flow at such boundaries. This can be achieved by partially incorporating the imaginary 
vertical distances corresponding to such boundaries through the use of certain coefficients 
in the segmental wetted perimeters. This can be achieved in many ways. This study 
investigates two different approaches. The first approach is based on Poesy’s (1967) 
approach. In tliis approach vertical imaginary boundaries are included to calculate the 
wetted perimeter in the middle sub section but not in the other sub sections. The second 
approach includes the imaginary vertical boundaries into all segments of a trapezoidal 
channel cross-section for the purpose of the calculation of the wetted perimeter. The 
calculations of the wetted perimeters according to the two approaches are described below. 

4.7.1 Approach-I: Model Q1 

This approach is based on Posey (1967) considerations for compound sections for 
calculating wetted perimeter. As per Posey (1967) criterion, the imaginary vertical 
boundaries are considered for the middle portion only and neglected completely in 
calculations of wetted perimeters for shallower portions. The incorporation of the 
imaginary vertical boundaries into the wetted perimeter calculations in the middle segment 
of the trapezoidal cross-section as per Posey (1967) criterion results in the following 
expressions for the segmental and total wetted perimeters. 
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(4.60) 

(4.61) 

(4.62) 

(4.63) 


where P,vi> Pw27 Pwi the wetted perimeters (m) in sub-sections 1, 2, 3, respectively; 
is the total wetted perimeter (m); y is the flow depth (m); zi and Z 2 are the side slopes for 
the right and left sides of the trapezoidal cross-section, respectively; and k\ and kz are the 
constants that account for the shear stresses at the vertical imaginary boundaries. 


4.7.2 Approach-II: Model Q2 


The Posey (1967) criterion is normally used for compoxmd channel sections consisting of a 
deeper middle portion and shallower portions on either side of the deeper section. Since 
the flow depth in the shallower portion in compound sections is small, the consideration of 
imaginary vertical boundary in wetted perimeter calculations for deeper segments only 
would not cause significant errors. However, in a composite trapezoidal channel cross- 
section, the maximum flow depth in all the segments is same. Therefore, neglecting the 
imaginary vertical boundary in wetted perimeter calculations for the right and left 
segments may be a source of error. A better approach may be to incorporate the imaginary 
vertical boundary into the wetted perimeter calculations for all segments of a composite 
trapezoidal section. This has been the motivation for the development of Approach-II in 
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this study. The incorporation of the imaginary vertical boundaries into the wetted 
perimeter calculations in all segments of the composite trapezoidal cross-section results in 
the following expressions for the segmental and total wetted perimeters. 


iy) + Ky 

(4.64) 


(4.65) 


(4.66) 

= ,J{zf+l)y,j(zl + l)+2(k,+k,) y + b 

(4.67) 


4.7.3 Determination of ky and kz 

The two approaches described above involve two unknown coefficients ki and kz that need 
to be determined. The question then arises is what the constants A-j and kz represent 
physically and how we can determine these. These constants should be functions of the 
velocity gradients since the velocity gradients at the boundaries are responsible for 
producing shear stresses and hence certain amount of resistance to flow in the vicinity of 
such imaginary boundaries. This can be achieved by considering the constants to be 
directly proportional to the differences in the velocity between right and middle segments, 
and left and middle segments, respectively. In other words, k\ can be assumed to be a 
function of (Vj-Vj) and kz can be assumed to be a function of (V3-V2). Please note that we 
should take absolute magnitude of the differences in the velocities in adjoining segments. 
Further, the maximum value of the constants, ^1 and kz, will be 1.0 corresponding to the 
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case when z\ zi 0.0 i.e. when the trapezoidal cross section reduces to a rectangular 
cross-section. Further, as the value of z, and/or increases, the value of h and kt is 
expected to decrease. This is because the effects of the rough side slopes on the flow field 
will tend to decrease as the side slopes go away from the imaginary vertical divides. The 
values of k\ and kt would attain their minimum values (say Arimm and fe min i when the side 
slopes become horizontal or the trapezoidal channel becomes rectangular without vertical 
boundaries. At such point of time the differences in velocities in the adjoining segments 
would dictate the values of these coefficients. Based on these concepts, the constants k\ 
and ki can be expressed by the following equations: 


Max{V,,V,) 


(4.68) 


Max{VM 


(4.69) 


Where V\ and Fz are the segmental velocities in the right and left segments corresponding 
to the side slopes of zi and zi, respectively, and Fs is the segmental velocity in the middle 
segment of the composite trapezoidal channel. It is to be noted that Fi, Fz, and F 3 will be 
varying depending upon the side slopes and roughness coefficients in different sub areas of 
the composite channel carrying given discharge. The constants k\ and ^2 would vary 
between Arimin and 1.0 and * 2 min and 1.0, respectively. The higher limiting condition (A:i = 
kt « 1 . 0 ) occurs when zi = zz = 0 . 0 , and Fi = Fz = 0 . 0 ; and the lower limiting condition {h 
* ^imin and kt = k%mx) occurs when zi = zz = <» (i.e. sides become horizontal). 
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I he formulation of the constants k\ and ki based on the principle of velocity gradients is 
presented by equations (4.68) and (4.69) above. However, the determination of ki and kj is 
not trivial since the segmental velocities Vi and V 2 on the right hand side of equations 
(4.68) and (4.69) will be functions of the wetted perimeters of various segments, which in 
turn depends upon the values of iti and ^2 themselves. Therefore, the determination of ki 
and ki would involve an iterative procedure. The determination of ki and k 2 through an 
iterative procedure is summarized in the following steps: 

1 . Determine the segmental velocities Vi, V 2 , and V 3 corresponding to the side slopes of 
Zi and Z 2 assuming the imaginary vertical boundaries as non-contributing to the 
wetted perimeter. 

2. Estimate the initial values of k\ and k 2 from equations (4.68) and (4.69), respectively. 

3. Using the estimated values of ki and k 2 , compute the revised wetted perimeters Pwu 
Pw 2 , and Pw 3 using equations (4.60) through (4.62) for Approach-I and equations 
(4.64) through (4.66) for Approach-II. 

4. Calculate the segmental areas (.^/’s) and determine the hydraulic radii (Ri = A\ I Pwd 
corresponding to the wetted perimeters computed in the previous step. 

5. Knowing the roughness coefficients («/, « 2 , ns) and the revised hydraulic radii (P,’s) 
calculated in the previous step, calculate the revised segmental velocities (Fi, V 2 , and 
Vi) using Manning’s equation. 

6. Repeat steps 2 through 5 until convergence is achieved in terms of the segmental 
velocities or total Q. 
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Oiicc the convergence is achieved, the final values of the constants, velocities, and wetted 
perimeters can be employed in the solution procedure through GA to determine the optimal 
dimensions of the composite channel as done earlier. The results from the optimization 
problem from model Q1 and model Q2 are presented in Table 5.13. 





CHAPTER 5 

Results and Discussion 

5.1 General 

Depending upon the method and parameters required to calculate the equivalent roughness, 
various models can be divided into three groups. The first group involves models A to D 
which requires the relative area of the subsections and hence depends on the method of 
subsection but not on the relative length of the wetted perimeter. These models do not 
represent physically reasonable situation especially if equivalent roughness (ric) is viewed 
from momentum perspective. The optimal design costs of composite channel in these 
models depend mainly upon the method of subsections. The second group involves 
models from E to H. They depend only upon the relative length of wetted perimeter and 
do not count the area of subsections. Hence, all models of this group have no effect on 
optimal design cost by method of subsection. Third group involves models from I to Q 
which depends on subsection areas and also on relative length of wetted perimeters. 
Hence, these models are expected to be more sensitive than others. 
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5.2 Results from Method-1 Scenario I 


Results liom Method- 1 Scenario-I that is division of total area into sub areas by vertical 
lines are given in Table 5.1. These models have been solved with GA. Results show that 
all the models of the first group are having approximately same cost varying from 
minimum cost 21.160 for Model C to maximum cost 21.241 for Model D. These costs are 
having much difference with the cost of models of another group. Hydraulic properties of 
all the models of this group are also similar. All the models of this group are having 
approximately same cross sectional area also varying from 28.925 m^ for Model A to 
29.092 for Model D. The models of the second group (Model E to Model H) are 
showing more or less same pattern. They have much difference with the models of another 
group but have quite similar values in the group ranging from minimum cost 22.702 for 
Model F to 23.026 for Model G. Cross sectional area of models of this group are also 
similar with minimum cross sectional area 31.454 for Model F to maximum area 31.893 
for Model G. Third group is more sensitive on the changes in the conditions as within this 
group any change in the ratio of sub area to total area relative length of wetted perimeter 
causes a lot of change in the search space hence the optimal design cost of the composite 
channel and seems to be more logical on momentum perspective. Hence, we can say this 
group is more diversified group giving optimal cost of composite channel ranging from 
14.379 for Model Q to 24.715 for Model I. These values are the best and the worst 
designing values of cost in all the groups. 
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To validate these results, these models have been also solved with the classical 
optimization technique, sequential quadratic programming (SQP) of MATLAB 7.01 
optimization tool box. The results are given in Table 5.4. This method is also giving 
similar values as in GA. Solutions by this method show that the cost values for the first 
group are ranging from 21.136 to 21.158 that is quite similar to the results with GA. For 
the second group these results range from minimum 22.685 for Model G to a maximum 
value 23.012 for Model H. Similar to the results of GA, this method is also giving the 
minimum cost of construction 14.341 for Model Q. 

Hence, it is clear from the results in Table 5.1 that Model Q is giving the minimum value 
of cost. I'his value is 30.96% less than the average value of all these models. Hence, we 
can say that trapezoidal channel designed by Model Q is the most efficient than the other 
models. Channel parameters (width, depth, and slopes of the two sides) and hydraulic 
parameters (area, wetted perimeter) are also given in Tables 5.1 for results with GA and in 
Table 5.4 for results with SQP. Results show that Model Q has minimum cross sectional 
area 1 8.723 m^ which is 34.27% less than the average value of all these models. 

5.3 Results from Method-1 Scenario II 

All the models have been solved for two other scenarios. Scenario II represent the 
designing of composite channel for all the models when there is present limited right of 
way in the field and we have to restrict the top width of channels. In this study, a constant 
value of top width equal to 7.5 m is taken as the upper limit to design the channel cross 
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section. For this purpose, an extra constraint represented by equation (4.55) was added to 
the optimization problem that restricts the top width to the maximum permissible top 
width. Results Irom Method -1 Scenario II that is models solved with method of division 
of total aica into sub areas by vertical lines and having an additional constraint, equation 
(4.55), solution by GA are given in Table 5.2. Results show that there is no change in the 
optimal design of the models of the first group as they all are having top width less than 
maximum permissible. They are having the same results as in Scenario I. Models of the 
second group are having higher values of the cost than in Scenario I but they have 
approximately similar values for all models ranging firom minimum value 22.726 for 
Model F to maximum value 23.066 for Model G. In this scenario also. Model Q has 
lowest value of cost in all the models that is 14.379 as in Scenario I. In this scenario, the 
average cost of all the models is higher than Scenario I and Model Q is 31.58% lesser than 
the average of all the models. In this scenario also. Model Q has least area of cross section 
(18.723 m^) that is 34.84% lesser than the average value of all these models. The results 
from classical method are given in Table 5.4. These results are approximately similar to 
the results of GA. 

5.4 Results from Method-1 Scenario HI 


All the models have been solved for another scenario that is restricted side slope. In this 
scenario optimization problem is subjected to some additional constraints, equations (4.56) 
to (4.59). Results solving all the models with GA are given in Table 5.3. In this method 
also, the models of first group are showing similar pattern as in the previous two scenarios. 
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These models have approximately same values of design cost having minimum value 
23.452 for Model C and maximum value 23.787 for Model B. In this scenario all the 
models ol' this group having higher cost than other two scenarios. Cross sectional area for 
these models are also approximately same. Models of second group are also showing 
similar pattern. They have design cost ranging from 24.330 for Model F to 24.615 for 
Model G. In this scenario, Model Q and Model L have least cost in all models that is 
20.931. This value is 1 1.20% lesser than the average of all the models. In this scenario, 
almost all ol the models have higher value of design cost than the two scenarios excluding 
Model I that has higher cost in Scenario 2. The results from all the models for Scenario III 
solved with classical method SQP are given in Table 5.6. These results are also giving 
approximately same values as calculated with the help of GA. 

5.5 Results from Method-2 Scenario I 

Results from Method-2 Scenario I that is division of total area into sub areas by bisecting 
lines at every change in geometry or roughness are given in Table 5.7. These are the 
results with GA. In this method also, the models from the first group having 
approximately same value with minimum cost of 22.674 for Model C up to the maximum 
value 22.976 for Model B. Similarly, models of second group are also showing similar 
pattern having minimum value of cost 22.702 for Model F and maximum cost of 23.026 
for Model G. In this method also, Model Q is having the lowest cost that is 22.453 which 
is 2.66% lower than the average value of all the models. While comparing these results 
with that of Method- 1 Scenario I of Table 5.1 we are getting some very interesting results. 
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Models of first group are giving slightly higher cost that means this method of division of 
total area into sub areas is giving over designed values for the models of this group. There 
is no change in the cost of models of second group. It is due to the fact that they are not 
al looted by the method of sub area division as explained earlier. For the models of the 
third group, this method is giving very high value of cost for all the models. One another 
conclusion that can be made while comparing these results with that of Method-1 is that 
excluding Model I (which is giving very high value in both the results) all the models are 
giving approximately similar values in Method-2 in comparison to Method- 1 as standard 
deviation lor cost of all the models for this method is 1.056 as that for all the models in 
Method- 1 was 2.964. Mence, we can conclude that the results for all models in this method 
are near to the average value as compared to Method- 1. For the Model Q which is 
perfonning in the best manner in the two methods, in the second method cross sectional 
area is 30.963 m^ which is 1.654 times higher than that in Method-1, the results all the 
models for Method-2 Scenario I solved with classical method (SQP are given in Table 
5.10. While comparing with GA results, these results are giving approximately similar 
values. 

5.6 Results from Method-2 Scenario II 

Results for Method-2 Scenario 11 are given in Table 5.8. While comparing to Method-2 
Scenario 1, Table 5.7, all the models of the first, second and third groups are having 
slightly higher values of cost as compared to Scenario I. Model I has the highest value of 
cost 27. 139 as Model Q has lowest one 22.885. When these results are compared with that 
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of Method- 1 Scenario 2, Table 5.2 the following conclusions can be made. Models of 
group one have slightly higher values of cross sectional area and also the design cost. 
There is no change in the results of the second group. Members of the third group are 
giving very high cost. These results from all the models for Method-2 Scenario II solved 
with classical method are given in Table 5.11. These results are approximately similar to 
that of GA in all respects. 

5.7 Results from Method>2 Scenario III 

Results for all the models for Method-2 Scenario III are given in the Table 5.9. These are 
the results with GA. All the models are having approximately similar cost. Least cost 
Model is Model Q with a cost of 23.524 and highest cost Model is Model I with a cost of 
29.076. While comparing within the results of the Table 5.9, Model Q has 4.22% lesser 
than the average value of all the models. While comparing with the other two scenarios of 
Method-2, all the models are having higher value than other two scenarios. When these 
results are compared with the results of Method-1 Scenario III, Table 5.3 the following 
interesting result may be noted. Costs of models of the first groups in this method are 
slightly higher than that of Method- 1. Cost of second group of models remains unchanged 
and cost becomes very high in the case of models of third group. Method-2 scenario 3 for 
all the models have been solved with the classical method (SQP from MATLAB 7.01 
optimization toolbox) and the results are presented in Table 5.12. While comparing with 
the results of GA, all the models are having approximately similar values as in GA and 
hydraulic parameters are also similar. 
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5.8 Results and discussion for Model Q1 and Model Q2 


From the previous results, it is clear that the Model Q is the best Model for designing least 
cost composite channel cross section. Model Q represent the flow conditions in the 
channel more accurately due to following reasons: 

(a) Allowing spatial variation in segmental roughness coefficient as using a single 
equivalent roughness coefficient. 

(b) Allowing spatial variation in velocity across the composite channel cross section 

(c) Allowing for incorporating perimeter along imaginary line. 

In the above Methods A to P, for designing optimal cross section of composite channel, it 
was considered that the dividing line is a zero shear line and there is no resistance along 
this line. But Model Q is based upon the spatial variation of velocity across the cross 
section and from the observations of velocity in the different subsections, it is clear that 
there is large differences in velocities in two sub areas divided by vertical lines. Hence, 
there will be a large velocity gradient along the vertical line that is liable to modify the 
resistance to the flow. To calculate the effect of the vertical imaginary line on the optimal 
design cost of the composite trapezoidal channel, two approaches have been developed in 
this study. Cost of all three scenarios with approach 1 (Model Ql) and approach 2 (Model 
Q2) are given in Table 5.13. In this case, it is clear that as we are considering the 
resistance along the vertical imaginary line of the Model Q, resistance to flow will be 
increased and hence a larger size of channel will be required to fulfill the need of the 
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required design discharge. In approach 2, since, the resistance to flow is more, slightly 
higher values of designed parameters are expected that increases the cost of construction. 
The values of channel parameters for the two approaches are presented in Table 5.13. 
Results show that Model Q1 has cost value 21.384 for Scenario I and Scenario II which is 
slightly lower than that of Model Q2 which has cost value 21.443 for Scenario I and 
Scenario II. In case of Scenario III Model Q1 has cost value 24.010 that is slightly lower 
than cost value of Model Q2 for Scenario III which is 24.309. Since in the total cost 
between the two approaches is not significant, and the second approach is more accurate, it 
is recommended for optimal design of composite channel. About the channel parameters 
it is clear from the results Table 5.13 that Model Q2 has larger cross sectional area in all 
three scenarios than that of Model Q1 but it has lower value of wetted perimeter than 
Model Q1 . These may be due to the fact that we have considered larger resistances to flow 
at the vertical imaginary lines in Model Q2 than in Model Ql. Hence, it requires larger 
area to carry designed discharge. 
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Table 5.6 Results from various models for vertical line division method, scenario III solved vi^ith classical method 
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Chapter 6 

Summary and Conclusions 

6.1 Summary 

This study presents the findings of an investigation of using different roughness coefficient 
formulation based on tlie different assumptions on the optimal design of composite 
channels. For formulation of models, seventeen different methods of calculation of 
equivalent roughness given by the past researchers of hydraulic engineering have been 
used. All these equations have their own limitations and assumptions. Using these 
equations in optimization problem, these limitations and assumptions modify the search 
space of the problem. Basic problem formulation is similar to the problem investigated by 
Das (2000) and Jain et al. (2004). These models have been solved for optimal design of 
trapezoidal composite channel using genetic algorithm (GA) and with classical 
optimization technique or sequential quadratic programming (SQP) from MATLAB 7.01 
optimization toolbox. Three methodologies have been applied to investigate the impact of 
roughness on the optimal design of composite channels. The first methodology is to 
calculate equivalent roughness by dividing the total area into sub areas using imaginary 
vertical lines from the places of change in geometry or roughness, and the second 
methodology is to calculate the equivalent roughness bisection and the third methodology 
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is based on calculation of wetted perimeter along the imaginary vertical line. For this 
purpose two approaches has been used. The first approach is based on Poesy’s (1967) 
consideration in which, wetted perimeter along vertical imaginary line was calculated for 
middle portion only and completely ignored in the others and second is based on the 
assumption that wetted perimeter along vertical imaginary line should be calculated for all 
subsections. For these approaches also the optimization problem has been solved for all 
the scenarios. For first two methodologies all the models have been solved for the optimal 
design of composite trapezoidal channel. For all the models and these two methodologies, 
two more scenarios have been investigated involving restrictions on top width and side 
slopes that may be warranted in the field, such as limited right of way, side slope stability 
criteria etc. 

In an attempt to search for global solutions GA was applied to 20 different initial random 
populations for all the problems in this study. Though the global optimal solution can 
never be guaranteed, it was found that most of the search corresponding to a different 
initial random population of solutions resulted in similar optimal solutions indicating that 
the global optimal solution was probably reached. It was found that the results firom GA 
for all the models and all the scenarios are comparable to those obtained by classical 


methods that also validate these results. 
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6.2 Conclusions 


In this study, it was found that method of calculation of equivalent roughness can impact 
channel design considerably. The results from different modes having higher values of 
cost may not be correct. Model Q performed in the best manner amongst all the models 

I 

for two methods and all the scenarios. Solutions by this model for Method- 1 had 30.96% 
less cost than average cost of all the models and 34.27% less cross sectional area than the 
average cross sectional area of all the models. In case of Method-2, almost all the models 
had comparable cost but Model had least cost that was 2.66% lower than the average value 
of all the models. Also, in this case cost of all the models except models E to model H 
increased that showed that the results from this method are inferior to that of Method- 1. 
In case of Scenario II and Scenario III also the Model Q had least cost and cross sectional 
area. In Method-1 Scenario II it had 31.58% lesser cost and 34.84% lesser area than the 
average value of the all models. Similarly, in Method- 1 Scenario III it had 1 1.20% lesser 
cost than the average cost of all the models. In case of Method-2, for Scenario II and 
Scenario III this model had least cost and cross sectional area in all models. In the 
proposed methodology Model Q1 and Model Q2 had comparable cost in all scenarios, but 
in Model Q2 wetted perimeter along vertical imaginary lines was calculated in all 
segments that is a better approach and is recommended for optimal design of composite 
channel. 

Each method is unique with unique assumptions. The method should be selected with 
proper care given to the site conditions in terms of availability of raw material in local 
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areas and appropriate assumption for site should be identified. Then the method most 
close to the assumption should be selected. 

6.3 Future scope 


There are some another methods of division of area into sub areas. These methods can 
also be investigated. This study can be applied to the design of composite channels of 
different shapes like parabolic, triangular, semi circular etc. This study can also be used in 
the field problems having a main channel and a flood plain those have different roughness. 
In solving the optimization problem, for its practicability velocity constraint may be added. 
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Calculation of area of sub sections in bisecting line method 



Figure A.l: Calculation of area of sub section 

Let ABCD is flow area of the trapezoidal section. Let line BP is perpendicular to line 
DC and line BN is the line that is bisecting the angle ZAfiC. LetZCBM = 6 , hence 
the angle ZABN will be (90-6/2) as line BN is bisecting the angle ZABC. From the 
above Figure in triangle A BCP 

tane = — ■ (A.1) 

Similarly in triangle A BPN 

tan (90-0/2) = - = cot 0/2 
z 

Hence, tan 6 / 2 = z (A.2) 


From trigonometry 


tan0 


2 tan 0 / 2 
1-tan ^ 0/2 


(A.3) 
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Putting values of these two in equation (3) 

J_ -- " 

z, 1 - z ^ 

2 2 = 1 - 2 - 

it is a quadratic equation. Solving it we get 
z = --z, ±/z7+l 

Area of the sub section ABCN will be sum of area of A5CP and ABPN. 
be 

y\-^iz y) y = ± -/zT+l 

Hence, the area of the subsection A; will be 



Similarly, the area of A 2 can be calculated. It will be 



(A.4) 

(A.5) 

Hence it will 

(A.6) 

(A.7) 

(A.8) 



